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Recent measurements of the anisotropy of the upper critical field Bc2 on MgB2 single crystals 
have shown a puzzling strong temperature dependence. Here, we present a calculation of the upper 
critical field based on a detailed modeling of bandstructure calculations that takes into account both 
the unusual Fermi surface topology and the two gap nature of the superconducting order parameter. 
Our results show that the strong temperature dependence of the Bc2 anisotropy can be understood 
as an interplay of the dominating gap on the o"-band, which possesses a small c-axis component 
of the Fermi velocity, with the induced superconductivity on the vr-band possessing a large c-axis 
component of the Fermi velocity. We provide analytic formulas for the anisotropy ratio at T = 
and T = Tc and quantitatively predict the distortion of the vortex lattice based on our calculations. 

PACS numbers: 74.20.-z, 74.25. Op, 74.70.Ad 



Our understanding of the physical properties of the re- 
cently discovered superconductivity in MgB2 has made 
rapid progress since its discovery 0. Its high criti- 
cal temperature Tc = 39 K can be understood as aris- 
ing from strong conventional electron-phonon coupling 
to a high frequency phonon mode 0, y, 01 ■ Its pair- 
ing symmetry seems to be of conventional s-wave type 
0^1^. However, in contrast to conventional supercon- 
ductors, a number of recent experiments indicate that 
there exist tw o ga ps of different size in this compound 
mill 13 in Hi- This possibility is supported by band 
structure calculations, which have shown that the Fermi 
surface of this compound consists of four bands: two tr- 
type two-dimensional cylindrical hole sheets and two tt- 
type three dimensional tubular networks 0, Q| in good 
overall agreement with recent de Haas-van Alphen ex- 
periments [1^ . Microscopic calculations of the supercon- 
ducting gap based on band structure calculations have 
shown recently that indeed one should expect a big su- 
perconducting gap living on the a bands and a smaller 
one, induced by interband pairing interaction, living on 
the TT bands |3l llq. Impurity scattering, which in con- 
ventional superconductors tends to average out strongly 
differing gap values, in this case becomes ineffective, be- 
cause the a and tt bands possess different symmetries, 
making interband impurity scattering much weaker than 
intraband impurity scattering |l7l| . 

Recent measurements of the upper critical field Bc2, 
particularly its anisotropy, on single crystal MgB2 have 
shown a puzzling strong temperature dependence of the 
anisotropy ratio Bc2/^c2 between the ab-plane and the 
c-axis upper critical field 0, 0, |2l|. In conventional 
systems this ratio rarely changes by more than 10 to 20 
percent as a function of temperature. In MgB2 changes 
by more than a factor of 2 have been observed. It has 
been shown that such a strong temperature dependence 
of the anisotropy ratio can be obtained within a strongly 



anisotropic single gap model, possessing a small gap in 
c-axis direction and a big gap in ab-plane direction [2ll |. 
However, the gap anisotropy would have to be a factor 
of 10, which is too big as compared with experimental 
values. In addition, this scenario would be inconsistent 
with penetration depth studies which clearly indicate the 
presence of a small gap within the ab-plane [H |2l| . 

Here, we want to present first calculations of the up- 
per critical field Bc2, which take into account the multi- 
band Fermi surface structure seriously. We show that the 
strong temperature dependence of the anisotropy of Bc2 
can be traced back to the influence of the two topologi- 
cally very different Fermi surface types. While the cylin- 
drical Fermi surface sheets are dominating the behavior 
of Bc2 at low temperatures, leading to a large anisotropy, 
at temperatures approaching Tc the tt bands due to their 
much larger c-axis Fermi velocity start to play a more im- 
portant role, strongly reducing the Bc2 anisotropy. Thus, 
the strong temperature dependence of the Bc2 anisotropy 
appears as a cross-over from a low temperature a band 
dominated regime to a higher temperature mixed a and 
TT band regime. 

In order to include the multi-band Fermi surface struc- 
ture in the calculation of the upper critical field, we 
start our investigation from the fully momentum depen- 
dent multi-band formulation of the quasiclassical (Eilen- 
berger) theory of the upper critical field For that 

purpose we have to solve the linearized multi-band gap 
equation in the presence of an external magnetic field, 
which reads 
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Here, /„ is the anomalous Eilenberger propagator on the 
Fermi surface sheet denoted by a. Aa is the gap function 
on that Fermi surface sheet and A"" is the pairing inter- 
action, which becomes a matrix in the band indices. The 
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FIG. 1: (color). Fermi surface topology used for the cal- 
culation of Bc2 in this work. The vr-band is modeled by a 
half-torus, the cr-band by a distorted cylinder. 



brackets (• ■ ■ )q' denote a Fermi surface average over mo- 
mentum k' of the Fermi surface sheet a' . In Eq. we 
have already assumed that the gaps are isotropic s-wave 
on the Fermi surfaces, as indicated by experiment, but 
may have different values on different sheets. At Bc2 the 
anomalous Eilenberger propagator has to be determined 
from the linearized Eilenberger equation (for ujn > 0) 
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Here, vp.a is the (momentum dependent) Fermi velocity 
on Fermi surface sheet a and A the vector potential due 
to the magnetic field B = V x A. Eqs. and Q 
constitute an eigenvalue problem for (r) . For a given 
temperature T the solution Aaif) which solves Eqs. Q 
and ||2J) for the highest value of B determines Bc2 ■ 

Usually this eigenvalue problem is solved by a Lan- 
dau level expansion of Aaif) above the Abrikosov ground 
state of the vortex lattice. It has been shown recently, 
however, that a variational ansatz for Aa{r) correspond- 
ing to a distorted Abrikosov lattice leads to much better 
results for strongly anisotropic systems jilf and we will 
adopt that method here. In this method the ansatz reads 
Aa{r) = Aaip\{r) with ipKx^y) = i/jAie~'^x,e^y). Here, 
ijjA is the usual Abrikosov groundstate and r is a varia- 
tional parameter describing the distortion of the vortex 
lattice. The undistorted lattice corresponds to r = and 
T has to be determined by maximizing Bc2- Introducing 
this ansatz into Eqs. ^ and Q and using standard op- 
erator techniques (see e.g. Refs. 11^113) we are 
led to the following eigenvalue problem for in band 
space: 
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Here, A+ is the highest eigenvalue of A"" , which deter- 
mines Tc- The function la is given by the expression 

Jq sinhit \ /a 

(4) 

Here, w_fi,q and VF2,a are the components of the 
Fermi velocity perpendicular to the magnetic field B 
on Fermi surface a. According to Eq. ||2J) the crite- 
rion for Bc2 is that the highest eigenvalue of the matrix 
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becomes 1. Apparently, 



la' "jfr ) 

aX T — Tc this is fulfilled for Bc2 = 0, because la — > 0. 

Eqs. Q and allow to determine the tempera- 
ture dependence and angular dependence of Bc2 from 
microscopic grounds. The material parameters we need 
for the solution Rre the Fermi velocities vpcn 

{k) and the 

coupling matrix A"" , which can be taken from band- 
structure calculations. In order to simplify the analy- 
sis we restrict ourselves to two relevant bands, because 
the two cr-bands and the two 7r-bands are very simi- 
lar 3J. The cr-band can be described to a good ap- 
proximation by a cylindrical Fermi surface with a small 
c-axis hopping parameter. The 7r-band can be mod- 
eled by a half-torus as shown in Fig. ^ (for compar- 
ison see the F-point centered Fermi surfaces obtained 
from bandstructure calculations in Fig. 1 in Ref. jT^). 
For the Fermi velocities on these two Fermi surfaces we 
thus write VF,-R{0,(f>) — v p {cos 9 cos (j), cos 9 sin (j), sin 9) 
and WF,tT(fcc, </>) = wfct (cos 0, sin 0, Cc sin c/cc)- Here, 
is the azimuthal angle within the a5-plane, 9 S [f : ^] 
the polar angle of the torus, k^ the c-axis component 
of the momentum, and c the lattice constant in c- 
direction. The dimensionless parameter ec describes the 
small c-axis dispersion of the cylinder. The Fermi sur- 
face averages over the cylinder and the torus are then 
given hy {■ ■ ^ ^ fj^^^ dk^ J^"" d(p ■ ■ ■ and 



. — - - - , where k is the ratio of 
the two radii of the torus. The parameters of the Fermi 
velocities can be found from bandstructure calculations. 
Taking the values given in Ref. I^^, we determine vp-n ^ 
8.2 • lO^m/s, VFa = 4.4 • lO^m/s, and Cc = 0.23. The ratio 
of the two radii of the torus can be estimated from the 
Fermi surface crossings in Ref. to be about n = 0.25. 

For a 2 X 2 matrix A"" the criterion that the biggest 
eigenvalue of Eq. (O becomes 1 leads to the equation 
(1 - ri)l„ +ril^ + \nt = -A± {l„ + \nt) {l„ +lnt). Here, 
t — T/Tc, A- is the smaller eigenvalue of A"" and 
rj ~ \^-\^ is a dimensionless parameter describing the 
intcrband coupling strength of the two bands {rj — cor- 
responds to no coupling, rj = 0.5 to maximum coupling). 
From bandstructure calculations in Ref. Q the following 
effective matrix elements can be obtained: X'^'^ = 0.959, 
X-^"" = 0.222, X""" = 0.163, and A^'' = 0.278. From these 



we find A 



+ 



1.008, A_ 



0.228, and t] = 0.064. As 



it turns out the remaining parameter Aj 



only 
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FIG. 2: Temperature dependence of the anisotropy ratio 
r — 8^2/ Bc2 for the two-band model described in the text 
and different interband coupling strengths 77. Solid circles are 
experimental results taken from Lyard et al.|20(|. 
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FIG. 3: Anisotropy ratio F = B^l /^c2 as a function of inter- 
band coupling strength i] for T = (dashed line) and T — Tc 
(solid line). The dotted line shows the approximation for 
T — given in Eq. (|7|l. At Tc the anisotropy ratio is much 
more sensitive to small interband coupling strengths rj. 



weakly affects the results as soon as A_ is sufficiently 
smaller than A+, as is the case here. 

We have calculated Bc2 numerically for the material 
parameters above optimizing the parameter r such that 
Bc2 is maximized. In Fig. |21 we show our result for 
the anisotropy ratio F = / ^ function of tem- 
perature for different values of the interband coupling 
strength rj. For rj = 0, when there is no coupling between 
the two bands, the temperature dependence of F is deter- 
mined by the cylindrical cr-band because of its stronger 
pairing interaction. Here, F changes only by 20 percent, 
as one expects for an isotropic single gap superconduc- 
tor. When 77 is increased, however, the temperature de- 
pendence of F becomes more pronounced. Our result for 
the parameters given above is shown as the dashed line. 
Note, that our calculation is parameter free, relying only 
on the parameters given by bandstructure calculations. 
For comparison, also the experimental data by Lyard et 
al. are shown (solid circles). As will become clear 
below, the most important parameters determining the 
temperature dependence of F are the interband coupling 
strength rj and the c-axis dispersion parameter ec. If 
we allow these two parameters to vary somewhat, we can 
obtain an excellent fit of the experimental data. The dot- 
ted line shows our result for 77 — 0.121 and ec = 0.182, 
where the other parameters have been kept constant. We 
mention here that this set of parameters also gives a cor- 
respondingly good fit of the temperature dependences of 
B^2 ^nd B^2 separately including the upward curvature 
of B^2 that has been noted in the experiments (2^ . This 
is just an immediate consequence of the strong tempera- 
ture dependence of the anisotropy ratio and B^2 varying 
linearly near Tc because of an absence of vortex lattice 
distortion in this field direction. 

Apparently, the anisotropy ratio F turns out to be 
more sensitive to interband coupling rj close to Tc than 



at T = 0. In order to get a better physical understand- 
ing for this behavior, we want to discuss some limits, in 
which analytical solutions for Bc2 can be obtained. At 
first we consider the limit T Tc- In this case Eqs. (O 
and I0J can be solved exactly and we obtain for the an- 
gular dependence of Bc2 as a function of the angle /3 the 
magnetic field makes with the c-axis: 



Bc2{P) 
Bc2{,P = 0) 



= (cos2/3 + Asin2/3)-i/2 



(5) 



where A = 2 



'71 — \/ 2 \ — : — r-tf — Here, (Vc^)a is the 
average of the squared c-axis component of the Fermi 
velocity over Fermi surface a, while (w^hQ)a is the cor- 
responding in-plane quantity. The distortion parameter 

When we introduce the 
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T is given by e ^"^ = 

numerical values for vftt, vpa and k, given above we find 
for the anisotropy ratio 
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For Cc = 0.23 this result is shown in Fig. 13 as a func- 
tion of rj (solid line). From Eq. (jnj we see that at 
?/ = the anisotropy is given by the c-axis dispersion 
of the cylindrical Fermi surface and actually diverges 
for ec ^ 0. (Note, that this divergence could not have 
been obtained from a Landau level expansion above the 
Abrikosov ground state and our variational ansatz above 
is crucial for this result). Once the interband coupling ry 
is increased, the anisotropy quickly reduces. This reduc- 

tion becomes sizeable already, when 77 ~ {v^^'" )" ~ iTeg' 
Thus, it is the small c-axis Fermi velocity of the cylindri- 
cal Fermi surface as compared to the one of the 7r-band 
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FIG. 4: (color). Vortex lattice structure for magnetic field in 
c-axis direction (left panel) and in ab-plane direction (right 
panel) at zero temperature calculated from the two-band 
model described in the text. For comparison, the area of 
the unit cell has been kept fixed. 



which leads to a high sensitivity of T near Tc and a small 
interband coupling is already sufhcient to make the in- 
fluence of the TT-band visible. 

In the limit T — > the integration in Eq. can 

be performed and we find la{T, ^r) — h^^^(^2^^fi 



i (In (e^^r^i^^ ^ g2r^|^^^^^^ where \nj = 0.577 is Eu- 
ler's constant. Since the pairing interaction is dominant 
in the ci-band, we can employ two approximations: at 
first we set A_ — 0. This corresponds to assuming that 
superconductivity in the 7r-band is completely induced 
by interband coupling. As a second approximation we 
assume that the vortex lattice distortion r is also dom- 
inated by the cr-band. Then we can find t by just opti- 
mizing it at 77 = and then use it for evaluation at 77 > 0. 
With these two approximations the logarithmic averages 
in la can be performed and we finally get 



0) 1.246 
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The parameter r is found to be t = | Inec, which shows 
that at low temperatures there is no simple relation any- 
more between t and the anisotropy ratio. Also we men- 
tion here that at low temperatures Eq. (|SJ) does not hold 
anymore. In Fig. |3lEq. Q is shown for ec = 0.23 as the 
dotted line. In order to demonstrate the quality of these 
approximations also the fully numerical result without 
these approximations is shown as the dashed line. Ap- 
parently, now the anisotropy F is much less sensitive to 
interband coupling and only varies on a scale given by 
77 ~ 1/ (0.482 — 2 Inec). The reason for this is that ec 
now only comes in logarithmically. 

Eq. 10 shows that at T = the divergence for Cc ^ 
appears at all 77 < 0.5, while at T = Tc this divergence 
only shows up at 77 = 0. Numerically we observe that 
at finite 77 and Cc = there is a certain temperature at 
which S^l 

diverges. This temperature becomes smaller 
when 77 is increased. In Fig. 0] we show the distortion 
of the vortex lattice at high magnetic field and low tem- 
perature that we expect from our calculation. When the 



magnetic field is directed along the c-axis of the crystal a 
regular Abrikosov lattice is expected as shown in Fig. 0] 
(left panel). However, when the field is directed within 
the a6-plane we expect a distortion of e"^ = = 0.48 
as shown in Fig. ^ (right panel) . This prediction can be 
checked by neutron scattering or STM tunneling [30| . 

To summarize, we have calculated the anisotropy of 
the upper critical field for the two band Fermi surface 
topology shown in Fig. ^ Using parameters from band- 
structure calculations for MgB2 we find a strong temper- 
ature dependence of the upper critical field in agreement 
with recent measurements on MgB2 single crystals. Fine 
tuning of the parameters can yield a very good fit of the 
experimental data. We observe that the small c-axis dis- 
persion of the (T-band leads to a high sensitivity of the 
anisotropy ratio on the interband coupling strength near 
Tc, but not at T = 0. This suggests that the interplay 
of these two quantities leads to the strong temperature 
dependence of the upper critical field in MgB2. 
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